<?php

// funciones sistema MM1
function MM1SistemaOcupado($miu, $lamda)
{
   $Pw = $lamda/$miu;
   return $Pw;
}
function MM1SistemaVacio($miu, $lamda)
{
    $Pw=  MM1SistemaOcupado($miu, $lamda);
    $Po= 1-$Pw;
    return $Po;
}
function MM1ProbNSistema($miu, $lamda, $n)
{
    $P=$lamda/$miu;
    $Po=  MM1SistemaVacio($miu, $lamda);
    $Pn= pow($P, $n)*$Po;
    return $Pn;
}
function MM1Lq($miu, $lamda)
{
    $Lq=  pow($lamda, 2)/($miu*($miu-$lamda));
    
    return $Lq;
}
function MM1L($miu, $lamda)
{
    $L=$lamda/($miu-$lamda);
    
    return $L;
}
function MM1Wq($miu, $lamda)
{
    $Lq=  MM1Lq($miu, $lamda);
    $Wq= $Lq/$lamda;
    
    return $Wq;
}
function MM1W($miu, $lamda)
{
    $L=  MM1L($miu, $lamda);
    $W= $L/$lamda;
    
    return $W;
}
//Funciones Sistema MMS
function MMSSistemaOcupado($miu, $lamda, $S, $Po)
{
    $P= $lamda/$miu;
    $Sf=  factorial($S);
    $Pw=((pow($P, $S)*$S*$miu)/($Sf*(($S*$miu)-$lamda)))*$Po;
    
    return $Pw;
}
function factorial ($S) {
 
   
    $resultado = 1;
 
    if ($S > 0) {
 
        // De esta manera estamos calculando siempre n * (n - 1)!
        $resultado = $S * factorial(($S - 1));
    }
 
    // Devolvemos el resultado acumulado
    return $resultado;
}
function MMSLq($miu, $lamda, $S, $Po)
{
    $P= $lamda/$miu;
    $Pw=  MMSSistemaOcupado($miu, $lamda, $S, $Po);
    $Lq=$Pw*($P/($S-$P));
    
    return $Lq;
}
function MMSL($miu, $lamda, $S, $Po)
{
    $P= $lamda/$miu;
    $Lq=  MMSLq($miu, $lamda, $S, $Po);
    $L= $Lq + $P;
    return $L;
}
function MMSWq($miu, $lamda, $S, $Po)
{
    
    $Lq=  MMSLq($miu, $lamda, $S, $Po);
    $Wq= $Lq/$lamda;
    return $Wq;
}
function MMSW($miu, $lamda, $S, $Po)
{
    
    $L=  MMSL($miu, $lamda, $S, $Po);
    $W= $L/$lamda;
    return $W;
}

//perdida de erlang
function Perlang($miu, $lamda, $n)
{
    $p=$lamda/$miu;
    $dividendo;
    for($i=0;$i<$n;$i++)
    {
        $K=  factorial($i);
        $dividendo=$dividendo+((pow($p, $i))/$K);
    }
    $nfac=  factorial($n);
    $perdida=($p/$nfac)/$dividendo;
    return $perdida;
}

?>
